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Abstract 

Although quantum tunneling between phase space tori occurs, it is 
suppressed in the semiclassical limit h \ for the Schrodinger equation 
of a particle in R d u nder the influence of a smooth periodic potential. 

In particular this implies that the distribution of quantum group veloci- 
ties near energy E converges to the distribution of the classical asymptotic 
velocities near E, up to a term of the order 0(1/VE). 



1 Introduction 

Consider abstractly a self-adjoint operator H on its domain D(H) in a Hilbert 
space H. Then for e > one may call a pair 

$,E) eD(H) xR , |$|| = 1 , \\(H-E)$\\<e 

an e-quasimode In particular, eigenfunctions ip with eigenvalues E are 0- 
quasimodes. 

The existence of an er-quasimode (ip,E) implies that the operator H has 
spectrum a(H) in [E — e, E + e] . In particular we are sure to find an eigenvalue 
E in an interval [E — /x, E + fi] for e < /x, if we know that the spectrum in that 
interval is purely discrete. 

*CPT-CIMRS, Luminy Case 907, F-13288 Marseille Cedex 9, France, e-mail: asch@cpt.univ- 
mrs.fr 

tlvlax-Planck-lnstitute for Mathematics in the Sciences, Inselstr. 22-26, D-04103 Leipzig, 
Germany, e-mail: knauf@mis.mpg.de 



1 



If we know in addition that E is the only such eigenvalue, then, after choosing 
an appropriate phase for its normalized eigenfunction ip, we have 

11^11 <-• (1) 

However, due to near-degeneracies of a(H) there may be no eigenfunction 
of H near ip: 

The standard example is that of the Schrodinger operator on the line with the 
double well potential V(q) := (q— l) 2 (q + 1) 2 . Then for energies E < 1 one may 
construct h°°-quas\ modes localized in one or the other well, whereas all eigen- 
functions of H h exhibit parity. Here the near-degeneracy of the eigenenergies, 
which is of order 0(exp(— c/h)), is connected with tunneling between the two 
components of the energy shell of the classical system (see, e.g. Lazutkin, [§]). 

So in that case phase space tunneling survives the semiclassical limit, and 
one cannot confine a particle forever in a well. As a physical consequence one 
may mention the NH 3 microwave radiation. 

For higher degrees of freedom d these energy shell components generalize to 
invariant Lagrangian tori in phase space. If such invariant tori exist and one has 
some control over the bicharacteristic flow in their vicinity, then it is possible to 
construct e-quasimodes of high accuracy (e = h N ) and thus to extract precise 
spectral informations in the semiclassical limit h \ 0, 

However, because of near-degeneracies in the spectrum, in general one cannot 
draw any conclusion concerning the semiclassical eigen functions. 

In recent years refined epitactic methods allowed to produce semiconductors with 
periodic superlattices. The electrons in these periodic potentials have a small 
effective value of h, leading to interesting effects (see QT^J] ) . In this context it 
is important to know to which extent one may model the electronic behavior 
classically. 

This motivates our study of Schrodinger operators 

h 2 

H h = A + V on H:=L 2 (R d ) 

whose potential V G C°°(IR d ,IR) is periodic w.r.t. a regular lattice £ C W 1 . 
We may consider V as a function V : T — > R on the <i-torus T := R. d /C. 
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The £-invariance and the Bloch theorem imply that H h conjugates unitarily 
to the direct integral of the operators 

H\k):=\(D + hk) 2 + V on L 2 (T) (k 6 T*), (2) 

acting on 




where C* is the dual lattice with Brillouin zone T* := R d /C* and D := —ihV is 
the momentum operator. It follows that the spectrum consists of bands. Up to 
measure zero sets due to degeneracies, the eigenvalues E%(k) are analytic in k, 
and are non-constant, (see, e.g., Thomas ||10|| , Wilcox ||13|| , and Reed and Simon 
0). Thus the group velocity hT x VkE^{k) vanishes at most on a set of measure 
zero. 

On the other hand the symmetry E%(—k) = E%(k) of the band functions 
implies in the non-degenerate case that the group velocity vanishes for k = 
and the other 2 d — 1 fixed points of k i— > —k on T*. 

To see how this vanishing of the group velocity is connected with phase space 
tunneling, we consider the simplest case of d = 1 dimension (for d = 2 see also 

0)- 

In that case the energy shell T, E := H 1 (E) C V of the Hamiltonian function 

H(p, q) := \p 2 + V{q) on the phase space V := T*T, 

consists for energies E > V max := max, V(q) of two components, corresponding 
to ballistic motion to the right resp. to the left. These components are permuted 
by the time reversal transformation (p,q) i— ► (—p,q) on V. 

As the eigenfunction ^(fc) can be chosen to be real for k — 0, it is semi- 
classically equally concentrated on both (one-dimensional) tori corresponding to 
the energy E = E*(Q). 

The vanishing group velocity is one manifestation of that fact. Thus for 
k — 0, arbitrarily small values of h and large times t the quantum evolution 
exp(— iH h (k)t/h) and the classical flow $* : V — » V generated by H behave 
very differently. 

However we argue that for general quasimomenta k in T* phase space tun- 
neling is exceptional in the limit h\0. 

More specifically, we conjectured in that the quantum distribution of group 
velocities converges in the semiclassical limit to the classical one, see Conjecture 
|5TT| below. 
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We proved this in for the extreme cases of potentials leading to ergodic 
motion, and for separable potentials (which are the only known examples of 
periodic potentials leading to integrable motion). 

Here we show a similar statement for arbitrary smooth potentials and large 
energies, where KAM tori are known to dominate the phase space volume. 

After presenting the strategy in Sect. 2, we adapt in Sect. 3 Lazutkin's results 
on KAM-quasimodes to the present situation of a family H h {k) of differential 
operators. Thm. |?[7] contains our main result. It states that for large energies E 
a proportion 1 — 0(l/yE) of the eigenfunctions is semiclassically concentrated 
near a KAM torus. 

This then leads to a corresponding statement ( Thm. [5^) for the semiclassical 
distribution of group velocities, in accordance with the above conjecture. 

In a final section, we try to abstract our strategy. We argue that a mere exis- 
tence proof for a full set of ^-quasimodes with localized asymptotic velocities 
could imply the conjectured classical limit of the distribution of group velocities. 

Acknowledgments. We thank Ruedi Seiler and SFB 288, TU Berlin, for hos- 
pitality, for which J. A. also thanks MPI in Leipzig. 

2 Heuristics 

Before we turn to formal statements and proofs, we shortly describe the main 
ideas, starting with the following observation. 

Two given quasimodes associated to different KAM tori give rise to different 
expectations of the sub-principal symbol hk ■ D of the operator H h {k) defined 
in (|2|). Thus they can be separated energetically by varying the quasimomentum 
k, and for typical k in the Brillouin zone T* one should not have too many 
near-degeneracies of energies. 

Of course we must consider scales in order to make this argument work. In 
d dimensions the mean spacing E^ +1 (k) — E^(k) between the eigenvalues of 
H h {k) near E > V m i n is of the order h d . Thus a priori one must consider in a 
fixed energy interval about h~ d quasimodes which may lead to a near-degeneracy 
with a given quasimode. For ^-quasimodes we need an energy separation of at 
least h N . So iV should be larger than d. 

Such high precision KAM quasimodes are constructed in the book || by 
Lazutkin (see also the article by Thomas and Wassell for related results) 
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We apply this method after some straightforward adaptation to our family (|^) 
of differential operators. 

An important input for that construction consists in the refinement of KAM 
theory presented in the paper [7]] by Poschel. Roughly speaking one uses that the 
deviation of the Hamiltonian function H from an integrable one vanishes faster 
than any power of the phase space distance to the KAM set. In particular we 
may apply perturbative semiclassical techniques in some ^"-neighborhood of the 
set of KAM tori. 

A final remark concerns the phase space complement of the KAM set. In 
general we do not have any information over individual eigenfunctions 

and eigenvalues concentrating semiclassically in that region. 

In particular we cannot hope to lift near-degeneracies between such eigenval- 
ues and the energies of the KAM-quasimodes by changing the quasi-momentum. 
Moreover, if a quasimode is involved in such a near-degeneracy, there need not 
be any eigenfunction ip^{k) near to that quasimode. 

However, we can apply a box counting principle. We know from KAM theory 
that for large energies E the complement of the KAM set is of relative measure 

o(i/Ve). 

Then a Weyl argument implies that up to an exceptional set of relative size 
0(l/yE) the eigenvalues E^(k) near E are well-approximated by KAM quasi- 
modes. 

In the semiclassical limit these /i^-quasimodes (ip,E) are typically energeti- 
cally separated in the sense that the associated intervals [E — h N , E + h N ] are 
disjoint. We have at least one eigenvalue E^{k) in each such interval. Thus 
only an exceptional set of relative proportion 0(1/VE) of these intervals may 
contain more than one eigenvalue. 

So for typical k 6 T* most E^(k) are not near-degenerate, and thus the 
corresponding eigenfunctions ip„(k) are well approximated by quasimodes ip. 

3 KAM Estimates and Quasimodes 

In order to apply KAM theory to H with energies in 

I:=[(1-5)E,(1 + 5)E] (3) 

near E > 0, we change coordinates. So consider the d x d matrix L := 
(£i, . . . , id)/ (2%) of a basis (£i,...,£d) for the configuration space lattice C, 
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set V(ip) := V(Lip), denote by V := T*T the phase space over the standard 
torus 

t := R d /(2iiZ) d 
and define, using the matrix M := (L*L) _1 , the Hamiltonian 

H £ :V^R , H e (J,<p):=\(J,MJ)+eV(<p). 
Then for the diffeomorphism 

Me:V^T , (p,g)^(J,^):= (L'p/^^g) 

we have 

E ■ H l/E oM E = H, 

and the flow $* generated by H £ (w.r.t. the standard symplectic structure on V) 
is conjugate to the original flow, up to a change of time scale: 

^ E t oM E = M E ° ( 5> t (teK). 

$* becomes fully integrable for perturbation parameter e — 0. Namely 

*o( J o, <Po) = (Jo, Vo + Vo(Jo)t) 
with the frequency vector 

^J) ~ f. (4) 
co> is of independent variation, i.e. the matrix 

du (J) 



d.J 



M is of rank d. 



So we are in a situation to apply KAM theory, see 0. For 7 > and 
r > d — 1 we consider the Diophantine sets 

fi 7 := [u eR d \Vk e Z d \ {0} :\u-k\> l\\k\\- T ] . (5) 

These are asymptotically of full measure as 7 \ 0. 
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Denote the interval of new energies by I := [1 — 8, 1 + 8]. For e = the 
phase space region V £ := C V is of the form 

V = J°° x T. 

By KAM for |e| small there exist a smooth generating function S £ on J°° x T 
and a Hamiltonian K £ independent of the angle variables, with the following 
properties. 

• The frequency vector 

LU £ :J°°^R d , J°° ^ \7K £ (J°°) 

is nondegenerate, and coincides for e = with (|4|). 

• On the Cantor set := (w e )"~ 1 (fl 7 ) of actions 

H e (J°° - d v S e (J°°, <p),<p) = K £ (J°°) ((J 00 , ip) e J™ x t) . 

• The symplectomorphism 

f £ :J°°xf^V , (J 00 ,^ 00 ) (J,p) 

generated by J°°<^ — S £ (J°°,(p) is near to the identity. 

• For 7 = the set /C e := T s (^ xT)n? £ of <l*-invariant KAM tori is 
of Liouville measure 

vol(£ £ ) >vol(P £ )- (l-O(Vi)). 

• The difference between the non-integrable Hamiltonian function H £ (x) and 
the integrable Hamiltonian K £ oT e _1 (x) vanishes faster than any power of 
the distance dist(x,/C £ ) from the invariant tori, and the same is true for 
any derivatives. 
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These statements imply corresponding results about the symplectic map T for 
the generating function S := \fESx/E o M. E and the integrable Hamiltonian 

K := E ■ Kx/eoMe 

J™ E := v^)- 1 ^ 
and the subset 

K T := Mi\£ s ) C Pj := Ho\l). (6) 
of KAM tori for the flow $*. In particular, 

voKXj) > vol(Pj) • (1 - OiE- 1 ' 2 )) . (7) 

Turning to quantum mechanics, the following theorem was essentially proven 
by Lazutkin in 0. 

Theorem 3.1 Letr > 2d in (ft), < H < 1 andk G T*. Define for a G (1, ^) 

Aftfc) := {£* G £* | dist(h(t + k),J™ E ) < H a }. 

Then for (3 := 1 - ad/(r - d) > 

2. (27r^) rf |A^(A;)| = l/C/l + O (HP). (8) 

Furthermore for N £ N, h small enough and t G A^(k) there exists a h N+l - 
quasimode (Ep(k), ^(k)). It follows that: 

2. there is an eigenvalue E h (k) of H h (k) with 

\E\k) - E%(k)\ < H N+1 ; 

3. for the spectral projection P on (E h (k) - h p , E\k) + h p ) it holds 

IKi-p^AOH <h N+1 - p . 

4. Let N > 2d + 2 and < p < N + 1 - d. Then Ve > 3a such that the 
dimension N of the space of all these quasimodes projected to the spectral 
subspace of\J et ,(E n (k) - K p ,E h (k) + K p ) meets the estimate 

(2Kh) d N=\1C I \ + 0(h l - £ ). (9) 
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Proof. This is essentially Theorem 41.10 in [§]. We specialize some formal 
aspects to our case - i.e. the configuration manifold has no boundary and the 
invariant Lagrangian tori are diffeomorphically projecting to the configuration 
torus, so that we do not need a Maslov operator, 
(ad 1): This is Lazutkin's Proposition 40.2. 

(ad 2 and 3): Let E be so large, that the KAM results hold true. The Ansatz 
for the quasimodes is: 

N+l N 

E\k) = h j Ej(k), $ h (k){q) = e US(k, q )-h(k, q )) ^ h j Aj(k, q) (10) 

3=0 j=0 

with 

#,)eC°°(T) , 

S(k, .) G C°°(]R d ) , S{k, q + £) - S{k, q) - h(k, t) E 2ttZ (£ G C). (11) 
Employing the operators 

T k ■= ~\{d q Sd q + d q d q S) = -i((V q S) ■ V, + \AS)), 
one computes 

N 

(\{d q SY + V-E«)YVA 3 + 

j=0 

N N+l 2N min(iVj) 

Y K +1 T k A 3 - \ Y h j+1 AA j _ 1 - V +1 Y E ^-^ ( 12 ) 

3=0 j=l j=0 Z=max(0,j-A r ) 

and is led to consider the equations 

H{d q S(k, q),q)- E {k) = O (h°°) (SC)-i 
and for < j < N, with A_i := 

T k A 3 (k) - \AA^ x {k) + ELo E J+ Uk)A(k) = O (h°°) (SCy 
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with the boundary conditions specified in fln|). 

The first step is to find a solution of the Hamilton-Jacobi equation [SC)-\. 
By KAM we know that there exists K G C 00 ^ 00 ), S G C 00 ^ 00 x T) such that 
not only 

H(P - OgSiP, q),q)= K(P) on (dpK)- 1 ^) x T 

but 

H(P - d q S(P, q),q)= K(P) + O (dist(P, (d P K)-\Q y )D (13) 

on J°° x T with all derivatives. Now set 

S(P,q) := Pq-S(P,q). (14) 
S{t,k,q) := S{h{t + k),q), E (t, k) := K(h(k + t)) 

then defines a solution of (SC)-i. 

Using the same strategy the transport equations (SC)j are now solved in two 
steps: first solve the corresponding equation indexed by P approximatively near 
a KAM torus, then replace P by h(t + k) for t G A^k) and exploit flatness of 
the functions. 

Let E be so large that d qP S(P,q) is non-degenerate. | det d qP S(P, q)\ dq is 
(the coordinate representation of) an invariant measure on a KAM torus P = 
const. So, with Tp denoting the transport operator with respect to S(P,q): 

(d g Sd g + (AS))\ det d 2 qP S(P, q)\ = 0^T P J | det d 2 qP S(P,q)\ = 0. 

V v ' 

For arbitrary P it follows that T P A (P,g) = O (dist(P, (dpiTp 1 ^)) 00 ) so 

A (q, t, k) := A (q, H(t + k)), E 1 := (15) 
satisfy (SC) for t G Aftife). 

By ( |I5| ) we may now suppose that A (P, q), E X {P) . . . Aj(P, q), E j+1 (P) meet 
(T P A f - lAA r ^ + ZCo 1 Ef+^A^P, q) = O (dist(P, (d P K)- l (n y ))°°) . 
Then the structure of the equation for Aj + i,Ej +2 is 

T P A j+1 (P, q) = f(P, q) + E j+2 {P)A {P, q). (16) 
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This is satisfied for P G (9piT) _1 (0 7 ) by 

E j+2 (P) :=- [ A -V(P,g(P,g)) dQ 

A j+1 (p, q (P,Q)):=MP, q (P,Q)) £ IH'm ^^ 

Here q(P,Q) is given by the canonical diffeomorphism T : (P,Q) i— > (p, g) 
generated by S(P,q), and g ^ g the Fourier- Transform 



£(f,P) := / g(P,Q)e-^^dQ. 
Jt 

Indeed, equation flip]) is equivalent to 

(-z0,S0,(V4b-i) = A o 7 + ^ +2 )(P ?) ^ 
-^o 1 ^! o **(0 9 S(P, g), g) f t=0 = (Ao 1 / + P i+2 )(P, g) 

where $* is the Hamiltonian flow of H. But Aq l A j+1 o $' o T _1 = A Mj+i o 
T- 1 o where **(P, Q) = (P, Q + d P Kt) is the flow generated by K. So 
equation ( |T5| ) is met by the above defined objects which are well defined and 
smooth if P labels a KAM torus and have a Whitney extension to J°° x T. So 
by the same argument as before 

A j+1 (t , k, q) := A j+1 (q, h(t + k)), E j+1 (t, k) := E J+1 (h(t + k)) 

satisfy (SC)j +1 . 

Define now with the functions so obtained the quasimode (xpf,(k), Ef*{k)) 
by the formula (|10]) with i>t*(k) normalized and the sum running up to iV; the 
sum for Ef*{k) runs up to iV + 1. We then have 

(H h (k) - E%{k))$%{k) = O (h N+2 ) 

so choosing h small enough we get the assertion. Items 2 and 3 follow by general 
considerations about quasimodes. 

(ad 4): To deduce (|9|) one has to estimate (^(k), ^^(k)), which is Lazutkin's 
Proposition 41.9. □ 

Remark 3.2 By it suffices to assume that the potential V G C Z (T, M) for 
/ G N large enough. 
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4 Approximation of Eigenfunctions 



Let the ft 27V -quasimodes {(ipf*{k), E$(k))} t *eAj(A;) be given by Thm. |3J] and 
denote by 

P*(k) (feA^)) 

the spectral projector for H h {k) and the interval [E%{k] - h N ,Ef*(k) + h N }. 
For each £* G Aj(&) there is a nearby eigenvalue 

with \E n n (k)-E%{k)\<h 2N . (17) 

So for h < ho we know in particular that dim(P^(A;)) > 1. 

But since the quasimode construction is only based on the KAM part of 
phase space, it does not suffice to know that the quasimode energies Ef*(k) are 
separated from each other to ensure that the eigenenergies are isolated. Thus 
we consider the subset 

TA^{k) := {t G QA^(k) | dim(P^(A;)) = 1} (k G TP). (18) 

of the index set 

QA^k) := jr G Aj(k) | \E%{k) - E$(k)\ > 2h N for g G A?(Jfe)\ {/ . } } , 
We obtain a map 

X k : QA^k) -> N 

by setting T k (t) := n for some n meeting (P7|). This map is one-to-one. 

Its restriction to PA^k) is uniquely defined, since for T G FAj(k) P^t(k) 
is the one-dimensional projector for the eigenfunction ipj k ^*)(k) of H h {k) whose 

eigenvalue E^^k) lies in - h N ,Ef,(k) + h N ). 

The index set QAj(k) of the separated quasimodes may be very small. For 
example it is even empty for k = in d = 1 dimensions, if h > is small enough. 
However, its mean cardinality 

over the Brillouin zone turns out to be asymptotic to 

(|^|)~(27r£r d vol(/C 7 ), 

with the KAM subset /Cj as defined in (|6|). This is the reason why indices in 
jFA^(fc) are abundant on the average; it holds: 
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Theorem 4.1 For t e JFAj(fc) , k e T* and a suitable choice of phase of the 
eigenfunction ipx k (t*)(k), 



(19) 



For N > d + 2 there is a (3 > such that for I := [(1 - 6)E, (1 + 5)E] with 
E > E th 

vol(/C 7 )-vol(/CJ)-C B (^) < (2irh) d < vol(/C 7 ) + C E (^),(20) 
with Kj :— Vi \ Ki. in particular 



{2irh) d (\FA?\) 
vol(Pj) 



(21) 



Remark 4.2 The Liouville measure of the thickened energy shell is of order 
vol(Pj) = c(5) • E d ' 2 • (1 + 0(l/E)). (22) 



Proof. Estimate (|TJ) follows from (|) and Def. (HI), since the (if)$(k), E$(k)) 
are h 2N -quas\ modes. 

The upper bound in (^) follows from the Lazutkin result (§) for |A/(A;)|. 

We claim that 



(27rfr) d (|£A^|) > vol(/Cj) - O(^). 
By (|8]) this follows from an estimate of the form 

(27ih) d (\A>}\gAi}\) = o(hP). 

But 

<|A/\SA/I)</ £ x(^W-^(A-) 

where x( x ) := 1 f° r \ x \ — ^ N ar| d otherwise. 
For E large and ft < h 

V - El(k))\ > \\V {E {^k) - E Q ^ 2) k))\ 



dk 

voir* 



(23) 
(24) 

(25) 



> W\£i - £ 2 \ > cte.H 2 



13 



uniformly for all k G T* and f i ^ l 2 6 Aj. Thus by the implicit function theorem 
the set of quasimomenta fceT* leading to a degeneracy 

El(k) = El(k) 

of quasi-energies forms a hypersurface, and 

lA^ (k) - mk) )^ = ° (hN ~ 2) - 

Since |AJ| is of order 0(h- d ), the r.h.s. of (||) is thus of order C(^ M+ ^ 2 ). 
So for iV > d + 2 + /3 estimate Q holds true, implying (|23]). 
We estimate the number 

(\^\) = (m\)-(m\fA ! }\) (26) 

from below by using ( |23|) and the relation 

\g^{k)\J=A^k)\ < \Ej(k) \X k (GAj(k))\ , (k G T) (27) 

where 

5*(*) := {n G N | Jj{(*) G /} 

is the index set of all eigenvalues in the interval /. 

Estimate ( p7|) follows by noting that its l.h.s. equals the number of intervals 

[E*(k)-h N ,E^(k) + h N ] for ieGA^k) 

containing two or more eigenvalues E%(k) (counted with multiplicity). By defi- 
nition of QAj(k) these intervals are disjoint, and we have 

El m {k) g \EUk) - n N M{k) + n% 

so that further eigenvalues must be indexed by an integer belonging to the set 
which appears on the r.h.s. of (|2"7|). 
The Weyl estimate 

{2nh) d |Sj(Jfe)| = vol(Vi) + 0{h) (k G T*) 

is uniform in k, since the slope of the band functions is bounded above by 

\V k E h n (k)\ < hy/2(E*(k) - V min ) 
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and thus of order H if E*(k) G / (see 0, Corr. 2.4). 
Thus the r.h.s. of (|7|) is bounded above by 

|S?(A;) \X fc (0\?(A:))| < (2vr^-' i vol(P / ) - |g/V?(fc)| - 0(ft 



Inserting that upper estimate for ( |27D in (|26| ) and using (^3|) proves the lower 
bound in 

Finally, estimate (|2"l| ) follows from (|2"fi| ) and the result 



vol(£J) 



(l/v^) , 



vol{Vi) 

see (0). □ 

5 Asymptotic Velocity 

As a consequence of Birkhoff's Ergodic Theorem for A-almost all xo E V 

1 f T 

v ± (x ) := lim — / p(t,x )dt 

T-^>±oo 1 J Q 

exist and are equal (A denoting the Liouville measure on V). In this case we 
set v(xq) := w ± (xo), and otherwise v(x ) := 0, thus defining the asymptotic 
velocity 

v : V -y M. d 

which is a A-measurable phase space function. 

We are particularly interested in the energy dependence of asymptotic velocity 
and thus introduce the energy-velocity map 

A := (H,v) : V — > R d+1 . (28) 

A is A-measurable and generates an image measure v := \A~ l on R d+1 . 

On the other hand (see [E|) for almost all fceT* the operator of asymptotic 
velocity 



(k) := lim - [ T e iHk W t (D + hk)e- iHh W t dt. 
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exists and is given by 

with the eigenprojections P^(k) of H h (k). 

The quantum asymptotic velocities are defined by 

f ^ _1 Vi.-E^(/c) , gradient exists 

K( k ) ■= { n *l ■ 

{ , otherwise. 

We equip the semiclassical phase space V h := N x T* with the semiclassical 
measure \ h := (27ch) d fii x /i 2 , where [L\ denotes counting measure on N and /12 
Haar measure on the Brillouin zone T*. 

In order to compare classical and quantum quantities, we introduce the 
energy-velocity map 

A h :V h ^ R d+1 with A h (n, k) := (E*(k), v*(k)) 

and the image measure v n := X h (A n )^ 1 . 

Example: For V = (free motion) v h = v independent of the value of h. 

In [§] we stated the following conjecture, which we proved for smooth V leading 
to integrable resp. to ergodic motion (see also for ergodic motions generated 
by Coulombic periodic V): 

Conjecture 5.1 For all C-periodic potentials V e C°°(R d , R) 

w* — lim v h = v 

(which means lim^\ L+i f(x)dv h (x) = J Rd +i f(x)du(x) for continuous func- 
tions f e C§(R d+1 , R) of compact support). 



Remark 5.2 One may also consider the stronger conjecture with continu- 
ous bounded test functions /, that is weak convergence in the language of 
probability theory. 

Here we obtain a statement which verifies the conjecture in the high energy limit. 
To this aim we introduce the ballistic scaling 

Me, v) : = E- d l 2 f{e/E, v/^E) (E > 0) 
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of a test function / e C§(M d+1 ,R), so that j\ = f. We notice that for V = 
we have u(E,v) = C ■ 5(E — \v 2 ) so that 



f E (x)du(x) = / f(x)du(x) (E > 0). 
The result is 

Theorem 5.3 For all f e C°(R d+1 ,R) with compact support in R+ x R d we 
have 



lim sup 

n\o 



f E {x)dv h (x) - / f E {x)dv(x) 



0(1/VE). (29) 



Proof. By our assumption on / there is an interval I of the form (|3|) so that 
I x M d strictly contains the support of f E . 

The index set of eigenenergies in I splits into the disjoint union 

H I R (fc)=Hf(A)UE 2 R (fc) with E^k):=X k (M^(k)). 
By the volume estimate (|22| ) and injectivity of 

(2nh) d (|E£|) = 0{E {d - 1)/2 ) + E {hP), 



so that 



(2ttW fE(E*(k),v*(k))dk = 0(l/VE) + E (hP). 

This leads to a contribution of order 0{\j\fE) to ( ^9|) , so that we need only 
estimate the contribution of E\. By (S) 

Um(2^(|A«\Af|> = for Af W : = {T <= A« W I «(r + *) 6 

So it suffices to consider the contribution of the index set 

E\ x {k) := T k {TK h j{k) n A*(fc)) c Eftfc). 
The result (^9|) then follows from the estimate 

u^(A;) = d P K(h(t + k)) + 0(h). (30) 
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for t G Tti}{k) n Af(Jfe) and n := J fc (f ) and the identity 

v(x) = dpK(P) (P ej°°,xe T({P} x T)) 
for the group velocity on the KAM tori which we both prove now. 

By definition ([18]) of Th!j(k), the eigenvalue E%(k) is non-degenerate so that 

*Z(k) = (il>*(k),v\k)rli(k)) (31) 
For (j) in the (^-invariant) domain of H h {k) and E G K. we have the estimate 

i||T? ft (AM| a < || - P)0|| II0H + \\V-E\\ ||0f (fc G T*). 

It follows from Theorem ( |4.1| ) that 

\\v\mi ie) (k) - = o (h n ) 

which implies for the expectation 

t(k) = (Vft*), (P + fi*)V5(*)> = + «0$t(A0) + O (/»"). (32) 

By construction of the quasimodes 

fa(k),(D + fUe)#£{k)} = J d q S(P,q) dfi P (q) + 0(h) (33) 

for P := h(t + k), S as defined in flTJ), and 

d 2 qP S(P,q) rfg 



rf/i P (g) := 



/ T ^, P S(P,<Z') dg'' 



Finally from the Hamilton-Jacobi equation, since the classical flow is ergodic on 
the invariant torus indexed by P, and since d[ip is the invariant measure in q 
coordinates, it holds 



v\x) 



[ d q S(P,q) dfi P (q) = dpK(P) (P ej°°,xe T({P} x T)). 

Thus (|30]) follows from Q, (f3|) and (0). □ 

Remark 5.4 Actually we have proven in addition to Theorem |3.1| that ip^ (k) 
lead to joint quasimodes of H h (k),v h (k), namely: 



\\{v h -d P K{h{t + k)))^Uk)\\ =0{h) (k G T*,t G jFAj(fc)). 
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6 Beyond KAM 



Theorem |5j3| gives a partial answer to Conjecture |5~|, based on the KAM region 



ICi C Vi. But what happens in the complement K,\ ? There the classical dynam- 
ics is very complicated in general, since one may encounter there further KAM 
tori (not predicted by the estimates), Cantori, elliptic and hyperbolic periodic 
orbits, large ergodic components etc. 

With the exception of the elliptic orbits, there is no direct generalization of the 
above KAM methods, and thus it seems hopeless to control the wavefunctions 
semiclassically supported in that region. However, as the following example 
shows, other methods may work. 

Example. Consider d = 2 dimensions. As shown in 0, in the presence of at 
least two geometrically distinct KAM tori the motion on E# is ballistic (v ^ 0). 
This is caused by the fact that these tori have codimension one in and thus 
confine the flow between them. We denote by IZi C Vj the phase space region 
enclosed by two nearby KAM tori (or rather families of such tori indexed by the 
energy in /). 

Using microlocal techniques, Shnirelman showed in [[3J the existence of a large 
number of quasimodes concentrated in TZj, see also 0. Now for large energy 
E the variation of v w.r.t. the restriction of Liouville measure to TZj is small 
in comparison with E. Thus by Egorov's Theorem the above quasimodes have 
group velocities near the classical v values (see also 0, Sect. 5). 

Different such regions IZi, however, have different classical asymptotic ve- 
locities. Thus one should be able to apply the heuristics developed in Sect. 2 to 
that case, too — without explicitly knowing the quasimodes. 

When trying to work on this kind of arguments, one is led to the paradoxical 
conclusion that sometimes it is more useful to know quasimodes (with certain 
additional properties) of an operator than to know its eigenfunctions. 
To explain this, consider the algebra generated by 

{H\k),$(k),...,vi(k)}, 

vf being the components of the operator of asymptotic velocity - which commute 
with H h (k) - and try to show the existence of joint quasimodes. Arguing along 
the lines of Sect. 2, such an existence proof could suffice to prove Conjecture [5J] 
in full generality. 
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